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FOREWORD 


This  report  consists  of  the  first  and  second  parts  of  the  three  part 

analytical  portion  of  the  High  Angle  of  Attack  Theory  (third  phase)  of  VRC's 

0 

program  of  developing  methods  for  assessing  the  non-uniform  flow  fields 
of  wing-propeller  slipstream  aerodynamics.  The  overall  program, 
consisting  of  four  (4)  phases  as  described  in  Ref.  1,  was  undertaken  to 
generalize  and  extend  the  basic  Rethorst  Lifting  Surface  Theory  (Ref.  2). 

The  first  phase  (Ref.  3)  was  comprised  of  an  application  of  the  basic 
Single  Jet  (Propeller  Slipstream)  Theory  (Ref.  2)  to  delineate  optimum 
finite  wing  planforma.  In  the  second  phase  (Ref.  1)  the  Single  Jet  Theory 
was  generalized  and  extended  into  a  Multiple  Jet  Theory.  The  generalizations 
and  refinements  of  the  second  phase  have  established  a  solid  foundation  for 
the  third  phase  (present  report). 

The  three  parts  of  the  third  phase  effort  are  comprised  of: 

1.  Wings  located  at  various  heights  in  the  jet 

2.  Highly  cambered  wings  as  used  in  deflected  slipstream 
V/STOL  arrangements 

3.  Tilt  wing  configurations  where  the  jet  is  at  an  angle  to  the 
free  stream  flow 

The  present  report  contains  the  basic  theorecticail  development  of 
the  first  two  parts  enumerated  above.  This  portion  of  the  third  phase  effort 
has  greatly  extended  the  applicability  of  the  analysis  by  encompassing 
deflected  slipstream  V/STOL  arrangements  currently  under  development. 


SUMMAJIY 


A  lifting  surface  theory  has  been  developed  for  wings  located 
at  arbitrary  heights  and  high  cuigles  of  attack  (up  to  the  inception  of 
flow  separation)  in  a  stream  containing  an  arbitrary  number  of  multiple 
jets  (propeller  slipstreams).  This  theory  extends  zind  generalizes  the 
formulation  of  T.  Y.  Wu  and  Richard  B.  Talmadge  (VRC  Report  No.  8) 
which  was  based  on  the  original  Rethorst  Lifting  Surface  Solution. 

The  present  theory  was  developed  by  first  analyzing  the  simpler 
single  jet  case  and  then  extending  the  analysis  to  encompass  multiple 
jets.  Thus,  the  analysis  is  systematically  presented  in  the  following 
order: 

1.  Wing  at  an  Arbitrary  Height  in  a  Single  Jet  -  a  method 
similar  to  that  employed  in  VRC  Report  No.  8  for  the  wing  located  along 
the  aixis  of  the  jet  was  used  to  extend  the  solution  to  a  wing  located  at 
any  height  in  the  jet. 

2.  Wing  at  a  High  Angle  of  Attack  in  a  Single  Jet  -  the  lifting 
surface  method  of  Weissinger  was  applied  to  chordwise  wing  sections  eaoh 
of  which  is  treated  in  accordance  with  its  height  in  the  jet  as  determined 
in  the  previous  step. 

3.  Wing  at  Arbitrary  Height  and  High  Angle  of  Attack  Extending 
through  Multiple  Jets  -  the  above  single  jet  analyses  were  extended  to  the 
case  of  a  wing  immersed  in  a  stream  containing  an  arbitrary  number  of  jets. 
The  jets  were  located  symmetrically  in  the  spanwise  direction. 

The  above  generalized  theory  encompasses  a  broad  spectrum  of 
V/STOL  aircraft.  The  prediction  of  wing-slipstream  interactions  afforded 
by  the  present  analysis  permits  optimization  of  such  V/STOL  aircraft. 
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I.  INTRODUCTION 


This  paper  provides  an  extension  to  the  work  of  T.  Y.  Wu  and 
Richard  B.  Talmadge,  who  considered  the  problem  of  a  wing  at  a 
small  angle  of  attack  placed  along  the  center  line  of  a  number  of  jets. 

Previous  attacks  on  this  problem  are  listed  and  discussed  by  Wu.  The 
effect  of  varying  the  height  of  a  wing  in  jet  slipstreams  has  not  been 
considered  previously  to  the  knowledge  of  the  author.  A  solution  is 
obtained  in  the  present  paper  for  a  wing  at  an  arbitrary  position  and  angle 
of  attack,  consistent  with  the  assumption  of  non- s eparated  flow,  in  a  stream 
containing  a  single  or  a  number  of  circular  jets  symmetrically  placed  along 
the  wing. 

The  problem  of  a  wing  at  an  arbitrary  height  in  a  single  jet  is  treated 
in  Section  II  by  similar  methods  to  those  used  by  Wu  for  the  wing  at  the  center 
of  the  jet.  In  this  latter  case  the  wing  is  represented  by  a  bound  vortex  line 
centered  along  the  quarter-chord  line  of  the  wing,  plus  a  trailing  vortex 
system,  the  boundary  condition  on  the  wing  being  satisfied  along  the  three- 
quarter-chord  line.  The  solution  for  a  wing  at  a  moderate  angle  of  attack  in 
a  single  jet  is  presented  in  Section  III  using  the  Weissinger  method  in  which 
the  wing  is  divided  chordwise  into  a  number  of  sections;  the  solution  for  each 
section  of  the  wing  is  then  given  by  the  relevant  solution  of  Section  II  corre¬ 
sponding  to  the  height  of  the  section.  In  this  way  a  bound  vortex  line  is 
placed  at  the  quarter-chord  point  of  each  section  and  the  boundary  conditions 
necessary  to  determine  the  vortex  strength  are  satisfied  at  the  three-quarte'r- 
chord  point  of  each  section.  The  problem  of  a  wing  at  moderate  inclinations 
extending  through  a  number  of  jets  is  attacked  in  Section  IV  by  extensions 
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from  the  single  jet  case  using  the  techniques  developed  by  Wu. 


The  solutions  described  above  are  derived  from  a  superposition  of 
infinitesimal  vortex  elements.  The  solution  for  each  vortex  element  is 
found  by  splitting  the  element  into  a  two-  and  a  three-dimensional  part.  The 
image  system  due  to  the  presence  of  a  circular  jet  in  the  main  stream  is 
easily  determined  for  the  two-dimensional  part.  A  similar  image  system 
for  the  three -dimen sionAi  part  is  assumed  plus  a  remainder  term,  which  is 
then  determined  by  the  boundary  conditions  on  the  jet  boundary.  A  feature  of 
this  approach  is  that  the  complicated  part  of  the  analysis  is  contained  in  the 
remainder  term,  which  for  the  examples  computed  by  Wu  so  far  has  contributed 
only  a  small  amount  to  the  total  lift.  No  justification  for  neglecting  this  part 
has  yet  been  found  however  and  it  remains  necessary  to  include  it  in  the 
present  treatment. 

Certain  simplifying  assiunptions  are  made  as  follows.  The  flow  is 
assumed  incompressible  and  inviscid,  in  which  case  a  velocity  potential  exists 
for  the  flow  inside  auid  outside  the  jets,  though  this  potential  need  not  be  con¬ 


tinuous  across  the  jet  boundary.  The  jet  boundary  may  thus  be  regarded  as  a 
circular  vortex  sheet.  The  slipstream  jets  are  assumed  to  be  circular  in 


cross  section  having  a  velocity  Vj  parallel  to  the  main  stream  of  velocity 


The  coefficient 


^  =  vy  V.  (1) 

is  formed.  The  jets  are  assumed  to  be  only  slightly  distorted  by  the  wing  so 
that  the  boundary  conditions  on  the  jet  boundary  may  be  applied  on  the  undis¬ 
turbed  circular  boundary.  These  boundary  conditions  are  the  kinematic 
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condition  that  the  flow  just  interior  and  exterior  to  the  jet  boundary  must  be 
parallel,  and  the  dynamic  condition  that  the  pressure  must  be  continuous 
across  the  jet  boundary.  Dealing  for  the  moment  with  a  single  jet  of  unit 
radius  whose  axis  is  taken  along  the  x-eixis,  the  velocities  outside  and  inside 
the  jet  are  taken  to  be 


V  =  V  i  +  V  ^ 
o  o  “o 


V.  =  V.  i  +  V 
J  J  J 


(2) 


where  (p^  and  <p^  denote  the  perturbation  velocity  potentials  outside  and  inside 
the  jet,  respectively.  The  boundary  conditions  on  the  jet  boundary  may  then  be 
expressed  as 


and 


on  q 


-I 


2^  2 
y  +  z 


=  1 


d(p.  dcp^ 


(3) 


(4) 


where  in  the  latter  a  linearized  Bernoulli  equation  has  been  used.  Eq.  4 
may  be  integrated  from  x  =  -  oo  to  give 

<Pj  =  on  r  =  1  (5) 

•I 

where  it  has  been  assumed  that  the  potentials  vemish  at  x  =  -  oo. 
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II.  SOLUTION  FOR  A  WING  AT  AN  ARBITRARY  HEIGHT 


IN  A  SINGLE  JET 


A.  FORMULATION 

The  analysis  presented  in  this  section  is  concerned  with  a  wing  repre¬ 
sented  by  a  bound  vortex  lying  along  the  line  given  by  x  =  0  and  z  =  -h  from 
y  =  -bj  to  y  =  b^,  and  its  trailing  vortex  system.  A  superposition  of  these 
solutions  along  the  chord  will  then  be  used  in  Section  3  to  represent  a  wing  at  a 
moderate  angle  of  attack  in  the  jet.  The  jet  is  located  with  its  axis  along  the 
x-axis  and  is  taken  to  have  unit  radius  (see  Fig.  1).  A  point  in  the  y-z  plane 
is  expressed  as 


r 


y  +  iz  =  qe 


ie 


(6) 


and  a  point  on  the  wing,  y  =  q  ,  is  given  by 

=  T|  -  ih  =  se”'°^  for  “h^  £  q  <  b^  (7) 

The  solution  for  an  infinitesimal  element  of  bound  vorticity  has  been 
given  by  von  Karman  (a  derivation  is  given  by  Wu).  The  result  for  a  bound 
vortex  element  of  strength  r(q)6q,  located  along  the  wing  x  =  0,  z  =  -h, 
for  q  <  y  <  q  +  6q  is 


b<p  = 


(8) 


where 
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6(^2=  -^r(Ti)6TiFj(y-Ti,z+h)  =  ■^r(T))6Ti 


(y-Ti)‘=^-K2^)‘ 


6</?3  =  r(ti)6tiF2(x  ,y-Ti ,  z+h) 


=  ^r(T,)6Ti 


{y-T|)'^+<z+h)  [x  -Ky-Ti)^+(z+ti)‘^]' 


The  velocity  potential  in  8  has  been  split  into  two  parts  6^  and  6^^,  repre¬ 
senting  the  two-  and  three-dimensional  parts  of  the  solution.  The  effect  of 
the  jet  on  the  solution  8  for  the  vortex  element  wiU  be  derived  in  parts  B,  C, 
and  D  of  this  section,  for  the  two^  and  three-dimensional  parts  of  the  solution, 
respectively.  The  solution  for  the  whole  wing  will  be  given  in  part  C  by  inte¬ 
grating  the  vorticity  distribution  r(T))  from  y=  -  bj^  to  y  =  b2. 


B.  IMAGE  SYSTEM  FOR  THE  TWO-DIMENSIONAL  PART  OF  THE  SOLUTION 


The  effect  of  the  jet,  represented  by  the  boundary  conditions  3  and  5, 

on  the  two-dimensional  part  6^2  ^  vortex  element  8  is  determined 

in  this  section.  A  simple  method  to  do  this  is  to  express  ^  e*  the 

difference  of  logarithmic  terms  representing  simple  line  vortices  situated 

very  close  together.  The  effect  of  the  jet  is  found  for  each  line  vortex  by 

choice  of  a  sxiitable  image  system  to  satisfy  the  boundary  conditions  3  and  5 

on  r  =  1,  the  boundary  of  the  jet.  Two  separate  cases  of  the  vortex  element 

2  2  2 

lying  outside  or  inside  the  jet  are  to  be  considered,  that  is  s  =  T]-i-h  >1 
or  s^  <  1,  respectively.  The  two-dimensional  part  of  the  solution  is  then 


-  5  - 


(i)  s  >  1.  i.  e.  for  vortex  element  outside  jet 


r(Ti)6Ti 


{fjCv-t),  z+h)  +  ^  ,z+-^  ) 


6  2  =  I  2r|hF ^(y  -  ,  z  +  •1}  for  q  >1 

r(q)6Ti  I  (l-€  2)Fj{y-Ti,  z+h)  -  nijFj(q,h)^  for  q  <  1 


(ii)  8  <  1,  i.  e.  for  vortex  element  inside  jet 


r(q)dii  I  (l-«  2)Fj(y-q,  z+h)|  for  q  >  1 

6<p^  = 

■—  r(q)dq  I  Fj(y-q,  z+h) - ^  ("(ri^-h^jF^Cy-  ,  z  +  ^  ) 

s'"  s  s 

+  2qhF2(y--^  ,  Z  +  ^  )  j  -  c  j^Fj(q,  h)|^  for  q  <  1 


where 


=  ^ 

1  7T^  ’ 

l+JX 


F3(y,z)  = 

y  +z 


(1-W‘ 


It  may  be  noted  that  the  terms  Fj^{q,  h)  in  the  solutions  above  are  constant 
terms  in  the  potential  and  hence  may  be  omitted;  their  inclusion  resulted  from 
using  the  boundary  condition  5  on  the  potential  rather  than  4  and  its  deriva¬ 


tive. 


The  two-dimensional  part  of  the  solution  corresponding  to  a  lifting 


vortex  line  of  given  circulation  distribution  r(y)  extending  from  ys  -  to 


-  6  - 


1 


I 

I 


y  =  may  now  be  obtained  from  9  and  10  by  linear  superposition.  Hence 
one  obtains 

b. 


4Tr  ^2  =  J  r(Ti)Fj(y-ii,  z+h)dTi  ”  «  2  J  z-rti)<iTi 


+  c 


z+  -^) 
s 


for  q  >  1  ,  (12) 


+  2T|hF o(y"  )1 

s  S'* 

^2  -h'  ^2  ) 

4Tr^2=  y  'r(q)Fj(y-q,  z+h)dq  "  «  2  ^  •^h'  | 

-  ‘1 1  [*’'  ^ 

V  _Vi  a  ^  a  a 


+  2qhF^(y-  IL^.  ,  z+-^  )J  dq 


for  q  <  1,  (13) 


where 


h'  =  ^1  - 


C.  THE  THREE-DIMENSIONAL  PART  OF  THE  SOLUTION  FOR  THE  VORTEX 
ELEMENT  LOCATED  OUTSIDE  THE  JET 

The  effect  of  the  jet  on  the  three-dimensional  part  5^^  vortex 

element  8  is  determined.  The  vortex  element  is  taken  to  be  located  outside 
the  jet,  i.  e.  s  >  1;  the  case  of  s  <  1  is  dealt  with  in  Section  Il-D.  The  image 
method  used  for  the  two-dimensional  part  cazmot  be  used  here  as  the  problem 
is  considerably  more  complicated.  The  technique- -introduced  by  Wu--of  using 
a  similar  image  system  as  the  two-dimensional  one,  plus  a  remainder  term 
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which  ia  then  to  be  determined  by  the  boundary  conditions  3  and  5,  will  be 
used.  Using  8  and  9,  the  proposed  form  for  the  vortex  element  is 


6^3=  r(Ti)6Ti  I  F2{x,  y-n,  z+h) 

+  ^  r(Ti^-h^)F2(x,y-  )+2rihF^(x,y-  ^  ,  z+-^  )1  | 

oo  a  a  ^  ' 


where 


6<p^  =  r(Ti)6Ti  (1-*  ^■*^)} 


^^(x,  y,z)=  ^  — 2 - 2^  21/2 


for  q  >  1  ,  (14) 

for  q  <  1  ,  (15) 


(16) 


and  Fp  has  been  defined  in  8.  Here  (p^  and  (pi  are  the  remainder  terms 

1 

left  in  the  potential  after  subtracting  the  F -functions  representing  a  convenient 
image  system:  (p'^,  <p'.  satisfy  Laplace's  equation  outside  and  inside  the  jet, 
respectively. 

Application  of  3  and  5  yields  the  following  boundary  conditions  on 
<Po  and  <p'. 

(<Pj-  = 

[(^^-h^)F2(x,y-  )+2qhF^(x,y-  ^  ,z  +  -^)j  | 


(17) 


( 


d(p'  d(p\ 


q  ■  ^  cJq 


-  *1  {  ^^(x.y-Ti.z+h) 

+ -^  r('n^-h^)F  ^(x,  y- ,  z  +  -^  )+  2qhF^(x,  y- ■3—  ,  z+-^  )!  I  (18) 
s'"  s  8  S  S  "*' 


-  8  - 


where 


y  =  q  cos  0,  z  =  q  sin  0 


(19) 


and  q  =  1  must  be  substituted  on  the  right-hand  side.  From  17  and  18,  or 
from  subsequent  forms  of  these  expressions,  it  may  be  shown  that  and 
<p.  and  their  derivatives  are  well-behaved  near  both  x  =  0  and  x  =  oo. 

Also  (p  and  w.  are  odd  functions  in  the  x- variable,  which  leads  to  the 
following  representations  of  these  harmonic  functions 


GO 


= 


II  {  sin  m0  j*  S<°>(k.s)K^(kq)  sin  kx  dk 
rn^  0  ^  ® 

+  cosm0r  (k,  s)K^(kq)8in  kxdk  I  for  q  >  1, 


(20) 


oo 


S^^(k,  s)I^(k,  q)8in  kx  dk 

m=  0 

+  cos  m0  ^  (k,  s)Im(l^q)sin  kx  dk|  for  q  <  1  . 


(21) 


The  functions  I  and  K  are  the  modified  Bessel  functions  of  the  first 
m  m 

and  second  kind,  respectively,  as  defined  in  Ref.  2.  is  used  in  the  solution 

•  » 

outside  the  jet,  and  is  used  inside  the  jet,  to  ensure  that  (p^  and  <p^  will 

be  regular  at  q  =  oo  and  q  =  0,  respectively.  The  super-indices  (o)  and 
(j)  chacterize  the  region  of  definition  of  the  functions  S  and  C,  i.  e.  outside 
or  inside  the  jet.  The  fiinctions  S  and  C  are  to  be  determined  from  the 
boundary  conditions  17  and  18. 

The  expressions  on  the  right-hand  sides  of  17  and  18  are  now  ex- 
panded  into  similar  series  and  integral  forms  used  in  representing  <p^  and 
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r 


(p!  above.  The  integral  representations  of  the  F-functions,  demonstrated  by 

F^{x.y.z)=  z  r  (£^y^+2^)-^/^d4,  (22) 

y+z  (x+y+z)'  ‘'0 


are  used  in  17  and  18  to  give  the  following  relations 


(<p'.-is(p'  )  =  pf  , /(sin  0  +h)  C  [4^+8^+1-29  008(0+0')) 

j  o  q-i  1  (  ^  ^ 

-  ^  hx  [  X  ^s  ^+ s  ^+ 1  -  28  cos(0+Q')]  ^  I 


(23) 


(d<p'o  / 

\aq  ■  -  ‘l{ 


x(  sin  ®+h)[  +  s  ^+ 1  -  28  cos  ( 0  +  a)]  ^ 


-  hx  [  s-cos(0+o')]  [  x^s  ^+  8  ^+  1  -  2s  cos  (0+a)]  ^ 

x^ 

(q  sin  0+ hq^)  [  4^+ q^9  ^+ 1-  2sq  cos  ( e+Q'))  ^d^l 


(24) 


in  which  q  =  1  has  to  be  substituted  on  the  right-hand  side. 

Tne  following  expansions  in  terms  of  Bessel  functions  are  used  in  23 

(4) 

and  24.  These  relations  are  given  in  Watson  or  may  be  deduced  from 
relations  there.  They  are 

x{x^+  s  ^+  1-  2s  cos  0)  —  r  sin  kx  K  (k  Ws  ^+ 1- 2s  cos  0)  dk  (  25) 

ir  J  0  o  ^ 

r  (|^+  s^+l-Zs  cos  9)  ^d4  =  —  C  ^  sin  kx  dk  ('  tK  (t^s^+l-2s  cos  0)  dt 

'''x  ’^''^0  '^Ik® 

(26) 
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xl  (lix)  K  f 
n  n 


'x)dx 


K  (wy)I 
n  ’ 


n  -1 


(jiy)  +  vK 


(^7) 


2. 


-  Zs  cos 


0) 


CXJ 

=  K  (ts)I  (t)  +  K  (ts)I  (t)  cos  m0 

o  o  '  m  m 

m=  1 

oo 

=  K  (t)I  (ts)  +  zY  K  (t)I  (ts)  cos  me 
o  o  m  m 

m=  1 


for  s  >1 

(28) 

for  s  <  i 


sin  0K  (t 

o 


Z 


+  1  -  2s  cos 


0) 


oo 

n=  1 
oo 

n=  1 


for  s  >  1 


(2'^) 


for  s  <  1 


Using  such  identities  as  above,  Eqs.  2i  and  24  may  now  be  written  in 


the  form 


(<p' -  yi(p'^) 


q=l 


IT 


sin  kx 


sin 


a[  K^(k)I^(|  )  -  K^(ks)I^(k)l 
s 


oo 

+  2  y  (  ^  r  K  (k)I  {  —  )  -  K  (ks)I  (k)l  cos  cc  sin  m(  oH  0) 
/  ,  Vks  '  m  m  s  m  m 

m=  1 


4  [  k'  (sk)I  (k)+  “Y  K  (k)I '  (  ^  ]  sin  a  cos  mlo+e)')  > 
m  m  c  m  ms  / ; 


(-^0) 
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dk 


=  — ^  sinkx|-k  sin  a  [l^(sk)Ij( s)  — ^  Kj^(k)Ij^(  ^  )  ] 

oo 

+  zY  (k[K'  (ks)I  (k)--LK  (k)I  (^)]3in  a  cos  m(a+e) 

/  ,  \m  m  iimms 

m=l  ® 

■  ^  f  cos  a  sin  m(cH-0)^|  dk 


I  I 

where  K  and  I  denote  the  derivative  of  K  and  I 

mm  mm 

The  boundary  conditions  3  and  5,  represented  by  30  and  31,  are 

now  in  a  form  suitable  for  comparison  with  similar  expressions  obtained  from 

the  expansions  20  cuid  21  assumed  for  the  solution.  The  S  and  C  coef- 

m  m 

ficients  in  20  and  21  may  then  be  determined.  Hence 


S^^Nk.s) 


,  2kl  r  K  (ks)  1 

j  \  ^  nri  Tin  nrx  I 

^  ®  l-€,k(I  k'  l'  )  J 

1mm  m  m 


cos  a  cos  ma 


■  “z  — : 

Ls 


2kl  l'  k'  (ks) 
m  m  m 


,k(I  K  +K  I  ) 
1mm  m  m 


d 


sin  a  sin  mQ'>  (3  2a) 


S^^^k,  s) 
m 


2«^(1-€  2)k 


I  k'  +K  i' 
mm  m  r 


l-€,k(I  K  +K  I  ) 
1mm  m  m 


(l^s) 

'  +K  l'  )  m 


cos  a  cos  ma 


+  K  (ks)  sin  a  sin  m( 
m 


(3  2b) 
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d°Nk,s)  = 

m 


.1S[- 


I  + 

m  3 


2kl  I  K  (ks) 
m  m 


iS.  IKS)  T 
n  m 

-  +K  l'  )  . 


cos  or  sm  mo 


l-«,k(I  K  +K  I  ) 
1mm  m  m 


r  1  ,  2kl  I  K  (ks)  T 

-  — f  I  (  — )  +  -  sin  a  cos  xna 

"  i-t,k(i  k'  +k  i’  )J 


,k(i  „  . 
1  mm 


m  m 


{32c) 


C^j^k.s)  = 

m 


2«^(1-€  pk 


I  k'  +  K  l' 
mm  — 


!-€  .k(I  K 
2  mm 


X-V.  i  / 

mm  I  ni  ^  V 

-  <  1 —  K  (ksjcos  a  sin  ma 

c’  +K  l’  )l  ks  m 


m  m 


K  (k8)sin  a  cos  mo  ? 
m'  ) 


(3  2d) 


Here,  the  argviment  for  the  Bessel  functions  of  argument  k  has  been  omitted 
for  brevity.  Eqa.  3  2  provide  the  coefficients  necessary  to  evaluate  ip^-  and 
(p\  from  20  and  21  and  hence  the  three-dimensional  part  of  the  solution  for 
the  vortex  element  lying  outside  the  jet  is  established. 


D.  THE  THREE-DIMENSIONAL  PART  OF  THE  SOLUTION  FOR  THE  VORTEX 
ELEMENT  LOCATED  INSIDE  THE  JET 

In  this  section  a  solution  is  developed  for  the  three-dimensional  vortex 
element  given  by  8  in  the  presence  of  a  circular  jet,  for  the  case  in  which  the 
element  is  located  inside  the  jet,  i.  e.  s  <  1.  Similar  techniques  to  U-C  are 
used.  The  vortex  element  is  represented  by 

5(^3=  r(r|)6r|  I  +  (l-€  y-'O,  z+ h)  |  for  q  >  1  (33) 

6^3=  -^r(r|)6Ti  I  +F2(x,y-Ti,  z+h)  -  [(q^-h^lF ^(x,  y-  ^  ) 

8  S  S 

+  2qhF^(x,y- ^  ^  for  q<l  (34) 

8  8  ^ 
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Here  <p^  and  (p\  are  the  remainder  terms  left  in  the  potential  after  sub¬ 
tracting  the  F-functions  representing  a  similar  image  system  to  the  two- 
dimensional  one  (see  10).  The  harmonic  functions  (p\  are  odd  in  x, 

and  regular  near  both  x  =  0  and  x  =  oo;  they  may  therefore  be  represented  by 
the  expansions 


oo  ^ 


m=0 


+  cos  m( 


i6  ^  D|°^(k,  s)K^(kq)sin  kx  dk|  for  q  >  1  (35) 

(p!  =  —  ^  isin  m0  r  T^'^^(k,  s)I  (kq)sin  kx  dk 
^  ^  J  Jo  m  *  '  m' 

10  ^  D^^(k,  s)I^(kq)sin  kx  dk|  for  q  <  1  (36) 


m=  0 


+  cos  m( 


The  potentials  33  and  34  may  be  substituted  in  the  boundary  condi¬ 
tions  3  and  5  to  give 

{(f  0  +  h)  r  [  ^  s  1-  2s  cos  (0+  O')]  ^  d| 

J  *1—  (  o  „ 


-  —  hx[x^s  s  1  - 2s  cos  ( 0-t- or)]  ^  (37) 


,  =  pe  I  x(  sin  0  +  h)  [x^+  s  1-  2s  cos  (  0-*  q)]  ^ 

'  ^  0=1 


-  hx  [  s  -  cos  ( 9+  O')]  f  x^s  s  1-  2s  cos  (  0-1-  a  )J  ^ 


x^ 

^  (q  sin  0  +  hq*^)  [  +  q  “^s  1-  23q  cos  (0+  a) 


1-3/2 

J 

(38) 
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in  which  q  =  1  has  to  be  substituted  on  the  right-hand  side.  Integral  repre¬ 
sentations  of  the  F-functions  similar  to  22  have  been  used  to  obtain  37  and 
38.  The  identities  25  -  29  are  used  to  express  37  and  38  in  the  expanded 
form 


2c .  oo  /  1  ic 

=— jp  J  sin  kx|sin  o[K^(k)I^(ks) - j  K^(  -  )  IJk)] 

0  S 

+  Z)  (Ik  (k)l'  (ks)  +  -^  k'  (  — )I  (k)l8in  ^cos  m(cH-0) 

yL  m  ^  m  s  rrr  J 


Fk  (  —  )I  (k)  -  K  {k)I  (ks)  1  cos  O'  sin  m( «+  0) )  f  dk 
ks  L  m  s  m  m  m  J  /  ^ 


(dtp'  9<^'\  i  r  1  W 

^  \  Sin  kx|-  k  Sin  o.[K^(k)I^(ks)  -^K^(|)I^(k)] 


oo 


(k)I  (ks) - -  K  (  — )I  (k)l  sin  a  cos  m(ol-  0) 

m  m  t  m  s  m  J 

s 


—  Fk  {  — )i'  (k)  +  K  (k)I  (ks)l  cos  a  sin  m(a+  0))  I  dk 
sLmsm  mmj  /; 


m  m 


The  boundary  conditions  39  and  40  are  now  used  to  determine  the 

I  I 

coefficients  T  and  D  of  the  expauisions  35  and  36  for  <p  amd  (D .  as 
mm  ^  o  ^  J 
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/V  I  K  +K_  I 

xj^^k.s)  =  2*j(1-€  2)k 


+  K  I 

mm  (  m  -  . 

—5 - r-  < 

K  +K  I  )  «  Ks  m 


cos  ma 


l-€,k(I  K  +K  I  ) 
1  mm  mm 


+  I  (ks)  sin  ti'sinmo' 
m 


} 


(41a) 


s) 

m 


=  Zc 


K  {-)+  - —  - T  COS  a 

1  '  ks  L  m  s  K  +K  I  J 


cos  ma 


1mm  m  m‘ 

,  2kK_K'  ll^(ks) 

K  (#)  - 


~~Z  m' s  ,  ,  , 
s  l-€ ,k( 


:  K  1  (ks)  I 
m  m  m  I 

(I  k'  +K  i'  )-• 


sin  a  sin  ma  /  (41b) 


ax  Tax  a 

mm  mm 


^(o), 


I  K  +K  I  , 

D-'(k.s)=  Z*  (l-«  )k  - 

l-€,k(I  K  +K  I  )  (  m 
1mm  mm 


s  a  sin  ma 


D^jJ(k.  3)  = 

m 


-  I  (ks)  sin  acos  ma/ 
m  } 

,  r  ,  2kK  k’  I  (ks)  1 

Im  K-.k.,  mmm 

liw  ^  cos  a  sin  ma 

11  ks  L  s  ^  K  +  K  I  )  -* 

1mm  m  m 


(41c) 


r  1  •  k  1  ) 

Ls^  ■"  =  l-,,k(I  K  +K  l'  )-l  > 

1mm  mm 


The  argument  for  the  Bessel  functions  of  argument  k  has  been  omitted.  The 
three-dimensional  part  of  the  solution  for  the  vortex  element  lying  inside  the 
jet  may  now  be  obtained  from  33,  34,  and  3  5,  36  together  with  equations  41. 
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E,  SOLUTION  FOR  A  LINE  VORTEX  DISTRIBUTION 


The  solutions  obtained  previously  for  an  arbitrarily  placed  vortex 
element  are  now  used  to  determine  the  solution  for  a  vortex  line  of  distri¬ 
bution  r{y).  The  vortex  line  lies  along  the  line  x  =  0,  z  =  -h ,  from 
y  =  -b  ^  to  y  =  b^.  The  two-  and  three-dimensional  parts  of  the  solution 
are  added  together,  using  the  relevant  solution  for  the  vortex  distribution 
inside  and  outside  the  jet,  to  give  for  the  perturbation  velocity  potential 
if?(x,  y ,  z) 


where 


4TT(p 


cp  =  (f^  for  q  >  1;  ip  =  (f ^  for  q  ^  1 


’2  ph' 

r(q)F(x,  y-q,  z+h)dr|  -  *  r(  q)F'(x  ,  y-q,  z+h)d  q 


-b 


1 


-h 


-  b 


+  € 


(q"^-h^)F(x,  y-  Z  +  ) 


+  ZqhG{x,  y-  ^  z  +  dq 


oo 


f 


oo 


1=  0 


K  (kq)sin  kx 
m 


(  (■/  r~^ 

s  sin  m0  1  \  '  i  / 


.(o). 


r(q)S''^'(k,  q)dq  +  T  q)T' "'(k  ,  qjdq 


.(o). 


-h 


+  cos  rn  0 


(  J  C/)r(q)d°>(k,q)dq  + 

'^-b  ,  ''h  '  -h 


r'(q)D  (k,  q)dii 


{-12a 


'  b. 


4TT^j  =  ^  r(Ti)F(x,  y-Ti,  z +h)d  n  -  w:  -i;  ,^r(Ti)F(x.  y-Ti,  z-fh)  dri 


€,  -^Tr'(^)r(n^-h^)F(x.y-  ^  .2+-^)dr, 

+  2TihG(x,y-  ,  z  +  — 2r)J  dr) 


GO  _ 

iY  r 

^  Jo 


I  ( kq) sin  kx 
m 


j,i„  „,e  [( r_  \  ') 

.  cos  -e[(5  "  +r  ‘) 


r(Ti)S^-»Nk.  r,)dTi  +  C  r(Ti)T'^Nk,  Ti)dTi 


.(j)  /u  4.  C 


+  \  1  r(  r))C  (k,  i)d  1  +  \  r(Ti)D  ■’  (k,  r))dr|  <  dk 

b  h'  '  *^-h'  "  J  ' 


(4^b) 


In  the  above  expressions,  the  functions  F  amd  G  are  given  by 
F(x,y,z)  =  F^(y,z)  +  F^(x,y,z)  =  ^  [l  +  x{x“+y‘^+z^) 


(43  a) 


y  +  z 


G(x,y,z)  =  F2{y,z)  +  F^{x,y,z)  =  ^ 

y  +  z 


[l  +  x(x^+y^+z^)"^/^]  (43b) 


respectively. 
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III.  SOLUTION  FOR  A  WING  AT  A  HIGH  ANGLE  OF  ATTACK 


IN  A  SINGLE  JET 

A.  FORMULATION 

A  solution  is  presented  here  which  determines  the  effect  of  a  single 
jet  slipstream  on  a  wing  held  at  an  angle  of  attack  which  is  not  small.  Angles 
of  attack  sufficient  to  cause  separation  of  the  flow  are  excluded.  This  section 
presents  the  solution  for  a  single  jet  located  centrally  along  the  span  of  a  wing 
whose  planform  is  taken  as  rectangular  and  whose  chord  is  untwisted.  This 
solution  will  be  extended  in  Section  IV  to  the  multiple  jet  case.  Further  ex¬ 
tensions  to  cases  of  non -rectangular  planform  and  small  twist  follow  from  the 
method  used  in  this  section  but  they  are  not  investigated. 

The  Weissinger  method  is  used  to  determine  the  flow  over  a  wing  at  a 
moderate  angle  of  incidence.  In  this  method,  the  wing  is  divided  chordwise  into 
a  number  of  sections.  The  flow  past  a  particular  section  is  approximated  by 
a  vortex  line  distribution  along  the  quarter-chord  line  of  the  section.  The 
solution  for  this  flow  is  obtained  from  Section  II,  taking  into  account  the 
position  in  the  jet  of  the  relevant  section.  The  total  solution  is  then  the  sum 
of  the  solutions  over  all  the  sections.  The  boundary  condition  of  zero  normal 
flow  at  the  wing  is  then  applied  at  the  three-quarter  chord  point  of  each  section 
and  at  a  sufficient  number  of  points  along  the  three-quarter  chord  line  to  com¬ 
pute  the  assumed  vorticity  distributions  to  reasonable  accuracy. 

The  wing  is  taken  with  its  leading  edge  along  the  line  x  =  0,  z  =  I, 
for  |y  I  <  b,  and  the  chord  line  (from  leading  to  trailing  edge)  at  an  angle  a 
to  the  stream,  e.  g.  see  Fig.  Z.  The  wing  shape  is  taken  to  have  the  form 

z'  =  S(x')  (44) 
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where  the  coordinates  x',  z'  run  from  the  leading  edge  along  the  chord  to 
the  trailing  edge  and  perpendicular  to  it,  respectively.  The  shape  of  the  wing 
in  the  x,  z  coordinates  is  taken  as 


z  =  H(x) 


(45) 


where  H(x)  .is  to  be  found  from  the  equations 

H(x)  =  f  -  x'  sin  a  +  S(x')  sin  o  , 
x  =  X  '  cos  a  +  S(x  ')  sin  a 

The  projection  of  the  chord  on  the  x-axis  is  given  by 
Cj^  =  c  cos  a 

where  c  is  the  chord  length  of  the  wing. 

The  wing  is  divided  into  M  sections.  Let 


A  section  is  defined  by  the  range 


(46) 


(47) 


(48) 


(n-l)d^  <  X  <  ndj^  for  n  =  1,  i,  .  .  .  ,  N  (49) 

A  vortex  line  distribution  is  placed  along  the  quarter-chord  line  of  each 
section,  that  is,  along  the  line  given  by 

X  =  (n-l)dj^  +  =  x^,  say,  (50) 

and 

z  =  H(x  )  =  -  h  say  .  (51) 

n  n  ’ 

The  boundary  condition  an  the  wing  will  be  evaluated  along  the  three-quarter 
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chord  line  of  each  section,  that  is,  along  the  line  given  by 


and 


X  =  (n-l)dj^j  + 


4  N  n 


say, 


(52) 


z  =  H(X^)  =  say  . 


(53) 


th 

The  vortex  line  distribution  at  the  quarter-chord  line  of  the  n 
section  is  represented  by  r(y)  =  r^(y).  The  solution  for  this  vortex  distri¬ 
bution  along  the  line  x  =  x^,  z  =  obtained  from  the  solution  given  in 

Section  II-E  by  substituting  (x-x^)  for  x,  h^  for  h,  and  rj^(y)  r(y). 

Call  this  solution  w  (x,y,z).  Hence 

^  n 


<P 


n 


<p(x-x  ,y,z;  r=r,h=h) 


(54) 


and  the  total  solution, 

N 


$ 


summing  over  all  the  sections  of  the  wing,  is  given  by 

(55) 


B.  FOURIER  EXPANSION  FOR  THE  CIRCULATION 

For  the  case  of  a  wing  square  in  planform  with  a  circular  jet  located 
along  the  mid-span,  it  is  convenient  to  adopt  the  following  notation  (see  Fig.  3). 
At  the  three-quarter  chord  line  of  each  section  of  the  wing  the  spam  is  divided 
into  two  regions:  R^  ^  outside  the  jet  and  R^  ^  inside  the  jet.  That  is, 

R,  :  I  y  I  <  h'  amd  R  :  h'  <  |y|<b  (56) 

l,n  I'l  n  o,n  nl^J 

where 

h'  =  Vl  -  hf  (57) 

n  »  n 
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The  angles  and  4^,  at  each  section  are  defined  by 

®^o,n  l,n  ’ 


y  =  b  cos  d* 


for  ^  ^  ^ 


I  I  I  ' 

y  =  b  cos  d*  =  h  cos  4*,  lor  y  <  h  (59) 

'  o,n  n  l,n  n 

The  circulation  functions  ^^e  expanded  as  Fourier  series  in 

terms  of  the  angles  id  »  'd  i  io  the  regions  R  and  R,  .  In  deciding 
®  ^o,n*  l,n  ®  o,n  l,n  ® 

on  the  particular  form  of  these  functions,  it  is  instructive  to  look  at  the  local 
lift  distributions.  By  the  Joukowsky  law,  the  lift  at  a  section  of  the  wing  is 
given  by 


!  (y)  =  f  V,  ,r  (y) 
n  ^  local  n  ^ 


The  lift  distribution  is  continuous  across  the  boundaries  of  the  jet.  Hence 


r  (h  -  ih  -  0)  =  pF  (h  -  ih  +  0) 
n  n  n  n  n  n 


r  (-h'  -  ih  +  0)  =  uF  (-h'  -  ih  -  0) 
n  n  n  *^000 


From  this  it  can  be  seen  that  the  circulation  function  F  (y)  has  a  disconti- 

n 

nuity  at  the  jet  boundary.  Also,  in  order  to  account  for  a  small  slipstream 
rotation,  odd  functions  of  y  in  F  are  allowed  inside  the  jet.  The  proposed 
expansion  for  rj^(y)  is 


F  (y)  =  4b V. 
n  J 


ZA 

y  sin 

Xd 

Lj 

o,n 

X=1 

F  +  A^^^ 

+  B 

0,  n  0, n 

0,  n 

in  R 


,  + 'V  sin  Xid  =  Ti  in  R, 

l,n  ^  X,n  o,n  l,n  l,n 
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in  which  A  and  A  may  depend  on  n;  that  is,  the  number  of  terms  taken  in  the 
expansions  may  vary  with  distance  of  the  section  downstream.  The  choice  of 
A  and  A'  will  be  determined  by  the  accuracy  required  in  the  solution  con¬ 
sistent  with  the  computing  procedure  adopted.  In  6Z  the  terms  and 

cos  4*  introduce  a  discontinuity  in  the  circulation  at  the  jet  boundary. 

The  slipstream  rotation  is  represented  by  the  term  cos  4^.  and  the 

^  ^  ’  o,n  l,n 

terms  involving 


gives 


Application  of  the  conditions  61  of  the  continuity  of  lift  at  the  jet  boundaries 


A-1 

=  (p.  -  1)  y  A^?^,  sin(2X+l)/3 

o,n  ^  /_j  ZX+l,n  ^n 

\=  0 


IX 

B^^^  =  (p  -  1)  y  A^.^  sin  ZXfS 

o,  n  ^  Z\,  n  ^  n 


where 


1  ^ 

/3  =  cos”  -ji 

^n  b 


On  neglect  of  slipstream  rotation  and  with  a  wing  symmetric  about  the  mid¬ 
span,  the  circulation  will  be  an  even  function  about  the  mid-span  which  indicates 

A^^[  =  A^°^  =  0  and  hence  =  0. 

Z\,n  2\,n  o,n 

The  total  lift  on  the  wing  is  given  by 

N 

L  =  y  L  (66) 

2^  n 
n=:  1 


where  is  the  lift  developed  by  each  section  of  the  wing.  That  is,  using 

60  and  62 
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=  4pbV 


.(I. 


VI  dy  +  ' 
o  o,n  ^  j 


v,r.  dyl 
J  1,  ^  f 


Substitution  of  the  expansions  6Z  for  F  and  F,  then  yields 

^  o ,  n  i,  n 


8pbV. 


h'(A<‘*  +iA">  )  +ibA'°' 
n  o.n  4  l.n  4  l.n 


A-1 

+  (p-l)bY  A^2\+l,n 
\=0 


sin  sin  2(X.+1)^ 


C.  DOWNWASH  DISTRIBUTION 

In  order  to  apply  the  boundary  condition  of  zero  normal  velocity  at  the 
wing,  it  is  necessary  to  calculate  the  downwash.  The  boundary  condition  will 
be  applied  along  the  three-quarter  chord  line  of  each  section  of  the  wing,  that 
is  along  x  =  X^,  z  =  (see  52  and  53),  The  downwash  velocity  at  x  = 


is  defined  as 


w{y)  = 


-[§?] 


x  =  X  ,z=H 
r  r 


r 

-  ^  t|P(X  -X  ,y,z;r=r  ,  h  =  h  )| 

0  7.  i_j  ^  r  n’-’’  n’  n 

J  z  =  H. 


where  (p  is  given  in  II-E.  The  downwash  is  split  into  its  even  and  odd 
functions  of  y,  which  are  denoted  by 


Ay)  = 


w  +  w  1 

o,ev  o,od 


w .  „  +  w .  , 

j,ev  j,od 


for  |y|  > 
for  |y|  <  h' 


1  -  H 
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The  subscripts  o  and  j  again  denote  the  regions  outside  and  inside  the  jet. 

Using  the  expansion  62  for  the  circulation  functions  in  the 

expression  69,  the  downwash  at  x  =  mav  be  evaluated  in  the  following 

form : 


-4Tr(w„  ^ 


o,ev^  ’^o.od'  "  Z-  I  "  3 

n=  1  ^ 


X  IksinGsininGK  (I'q)  *■ 
V  ni 


<  ,N  +  >N  +  }N  +  —  r  sin  k(X  -x  )  ^ 
jin  2n  3n  ^Jn  rn'/, 

P  m=  0 

(k  sin  G  sin  inG  K  (I'q)  —cos  6  cos  m6  K  (l<q)  ) 

\  ni  q  m  / 

I  I 

([l  b"-"  C  ‘ 

n  n 

+(  k  sin  6  cos  m6  K  (kq)-  —  cos  6  sin  m0  K  (kq)  1 
\  m  ^  q  m  ^  / 

t  t 

(  [r.b  "  "  I !■  ]  "  Ipnl"  >°  m  O'-  '> 


-  4'tT(  w  +  w . 

J  .  e  V  j  , 


N  /  CO  °° 

^)  =  )  I  ,N  +  +  .N  +  -  \  sin  k(X  -x  )  ) 

l,od  ZZ  I  1  4  n  5  n  ^Jn  ^  ^  Lj 

n=l  '  ^  m=0 


(71a) 


in  0  sin  m0  I  (kq)  +  — cos  0  cos  m0  I  (kq) 
m  ^  q  m  ^  , 


^k  sin  0  sin  m0  ! 

([€'•  d 


r  ,  r  (q)T|jNk.q)dq 


+  l  k  sin  0  cos  m0  I  (kq) - cos  0  sin  m0  I  (kq) 

\  m  q  m 


([C-d'-. 


(q)c[^^k,q)dr|  +  f  ,  F  J r, ) oj^^k.q  )dq^ 

(71b) 


-  25  - 


where  =  y^+  and  0  =  tan~^  have  to  be  substituted  on  the  right- 

hand  side.  The  N  functions  are  defined  by 


=  4b^v.y  -  .1;;  1 

1  '  X.,n  lo  2o  io 
X=1 

*  4bh;v.(A^;'„  [  +  5I  -  It  7I  +  e'  -  9'  ' 

lA  1. 

+  A<^>  [.I't  ,1",-  .l"  1) 

X,n  lo  2o  3o  / 


(7  2a) 


->  I  \  ^"3  '^~P‘ 

''j*  a  I  4  „  '  1^)3  "■*  2‘g  3‘3  "1 

-  4bh'  V.«  ,  [  .1  +  J  -  [-1  +  „I  -  -I  1 

nj2\o,n  4  5  6  o,n7  8  9 

2A'  V 

+  >  a1^)  [  +  ,1"  -  ,r]) 

K.n^lo  2o  3  o  / 


(72b) 


(  0  0 

/I1-2,,  V  .(o)  ,  n  .  .TT  .  .^n  .TT 

j  1  2L  1C  o  10  11  o  11  TT-p^ 

Xs  1 


/3  ^  » 

I  -  I  "  -  l"" 

12o  12^tt-3  19  o  19  IT-0  ) 

'^n  n  ) 


4ln  = 


wV  (  P  ft  ft 

.tb^v..,^  A<°'  ,I  "  +  «  3-  ,I„"  +  ,lZ„  -  ,1"  „ 

j  24^  X,n|lo  2  0  2  .-0^  3  0  3  ,-0^ 

\m1 


(72c) 


(72d) 
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2A 


cl  *  -4b^V.€,  )  A»  I 
5  n  J  1  Zu  X,  n  MO  3 


Y  .(«)  [ 

i  1 


X«1 


n-3  n-B  n-B 

+  I  ^'-I  *^-1 

11/3  12^3  19  ^  ^ 

n 


4bh'  V.€,  [,,I  +  ,  .1  -  1  +  [  wl  +  ,_1  -  'oI  -  o,ll 

n  j  1  \  o,n  13  14  15  20  o,n  16  17  18  21 


2A 

I 

X*1 


Ar'  f  +  111’'  -  1  -  iqI"  W 

X,n  10  o  11  o  12  o  19  o  / 


(1) 


(72e) 


In  the  above  expressions  a  dash  on  a  Bessel  function  denotes  differentiation 
with  respect  to  the  argument.  The  1  functions  are  defined  in  Appendix  1. 
The  even  and  odd  parts  in  the  N  functions  are  determined  by  taking  X  odd 
and  even,  respectively,  in  the  series  terms  together  with  the  terms  con¬ 
taining  in  the  even  part  and  in  the  odd  part.  The  remaining 

terms  on  the  right-hand  sides  of  7  2  may  be  split  into  their  even  and  odd 
parts  by  taking  the  even  part  to  contain  only  odd  values  of  m  in  the  series 
containing  and  T^  and  even  values  of  m  in  the  series  containing 

and  I^j^>  The  odd  part  is  then  given  by  the  remaining  terms  in  these  series. 


D.  BOUNDARY  CONDITIONS  ON  THE  WING 

A  slipstream  rotation  is  considered  which  has  the  following  properties. 
The  slipstream  may  have  a  rotational  velocity  component  about  its  own  axis. 
The  amgular  velocity  is  denoted  by  u>(q),  where  q  is  the  radial  distance  from 
the  axis;  u  is  taken  positive  if  the  rotation  is  counterclockwise  when  viewed 
from  the  rear.  It  is  assumed  that  ]  w  <<  V..  The  slipstream  rotation 

is  assumed  to  have  no  effect  outside  the  jet.  Its  effect  inside  the  jet  is  a 
change  in  the  fluid  direction  at  the  wing,  resulting  in  an  effective  variation  of 
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wing  incidence.  A  variation  of  the  angular  velocity  u(q)  in  the  downstream 
direction  may  be  included  using  the  analysis  presented  here,  but  this  will 
not  be  done. 

The  local  geometric  angle  of  attack  of  the  wing  along  the  three-quarter 
chord  line  of  the  r^^  section  of  the  wing  is  denoted  by  Q'y(y).  It  is  assumed 
that  the  geometric  wing  incidence  is  symmetric  in  y,  i.  e.  that  j.(y)  “ 

The  boundary  condition  of  the  tangency  of  the  fluid  velocity  at  the  wing  gives 

w(X^,y)  s  V^a^(Y)  for  |  y|  >  (73a) 

.y)  -  Vy^+*^r  |y|  ^  ^r 

r 

along  the  three-quarter  chord  line  of  each  section  of  the  wing,  that  is,  for 
r  =  1,.  .  .  ,N.  The  positive  square  root  is  taken  in  73b.  The  definition  70, 
splitting  the  downwash  w  into  its  even  and  odd  components  inside  and  outside 
the  jet,  may  be  used  to  achieve  a  similar  splitting  of  73  as 


for 

h' 

r 

<  y 

<  b 

(74a) 

o,od<^r'y)=  ° 

for 

1 

H 

r 

<  y 

<  b 

(74bf 

j,ev<^r-y)  = 

for 

0  < 

y  < 

1 

H 

r 

(74c) 

/y^+H^  )  sin 

JL 

^r 

for  0  <  y  <  H 
'  r 

(74d) 

Equations  71  and  74  provide  a  system  of  linear  algebraic  equations 

for  the  determination  of  the  coefficients  in 

X.,  n  X,  n  o,  n  o,  n 

the  Fourier  expansion  tZ  of  the  circulation.  The  boundary  condition  repre¬ 
sented  by  74  must  be  applied  at  sufficiently  many  points  y  along  the  wing 
to  obtain  a  determinable  system  of  equations  for  these  coefficients.  If  an 
equal  number  of  coefficients  is  used  at  each  section,  the  niimber  of  equations 
involved  will  be  N  times  that  of  the  single  section  wing  used  in  the  example 
treated  by  Wu.  A  feature  of  the  present  work,  however,  is  that  the  singularities 
present  in  the  Wu  solution  (where  the  boundary  condition  is  applied  on  z  =  0) 
do  not  appear.  This  is  due  to  the  presence  of  terms  of  the  type  (Hj.+ 
in  the  denominator  of  the  integrals  in  the  Appendix;  these  terms  never  vanish 
and  prevent  a  singularity  in  the  integral.  Under  these  circumstances  the 
constrains  on  the  coefficients  sufficient  to  exclude  such  singularities,  as 
adopted  by  Wu,  are  relaxed  here.  The  boundary  condition  74  must  be  applied 
at  a  compensating  number  of  points  y  along  the  wing  to  make  up  the  equations 
released  by  the  relaxed  constraints. 
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IV.  EXTENSION  TO  MULTIPLE  JETS 


A.  NOMENCLATURE 

A  system  of  2J  jet  slipstreams,  symmetrically  placed  with  respect 
to  the  mid  span,  is  considered.  The  wing  span  is  divided  into  (2J+1)  regions, 
denoted  by  R  ,  v  =  0,1,  .  .  .  ,  2J  which  are  defined  as  follows:  R  is  the 
region  outboard  of  all  the  jets,  1  region  inside  the  v  jet  pair, 

V  =  1,  2 . J,  counted  from  the  wing  tip  toward  the  center, 

region  between  the  v  and  the  (  v  +  1)^^  jet  pairs.  The  width  of  these  regions 
will  vary  in  the  downstream  direction  according  to  the  inclination  of  the  wing. 
The  boundaries  of  the  regions  along  the  y-axis  are  given  by  the  points 


b^=  b,b^,b^ . ‘^2J'^2J+1‘  " 


^-2J . ‘*^1’  '^o  ~ 


The  center  of  the  v  jet  is  along  the  line  defined  by  z  =  0  and 


y  =  +  a 

■'  —  V 


'=2v> 


The  wing  will  be  taken  to  have  a  square  plamform  and  to  have  no  span- 
wise  twist.  At  a  station  x  downstream,  therefore,  the  wing  will  be  a  distauice 

h  below  the  center-line  of  the  jets.  The  intersections  of  the  wing  with  the 

th  ' 

boundary  of  the  v  jet  at  the  station  x  is  given  by  z  =  -h  and  y  =  a  +  h  , 

'  I  1 

where  h  -  ^1-h  ,  which  we  denoted  by  y  =  c^^  and  y  =  ,  respectively. 

Hence  a  point  y  in  the  region  R^  may  be  found  from 

y  =  Cg  cos  4>g,  s=0,l,...,v  (75) 

At  the  point  y  =  c^ 

4.^  =  cos’M  .=  0,1 . V  (76) 

■  V  =  0,1,  .  .  .  ,  2J 
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Fig.  4  gives  a  guide  to  the  nomenclature  defined  above. 


B.  FOURIER  EXPANSION  FOR  THE  CIRCULATION 


The  Weissinger  approximation  will  be  used  to  determine  the  circulation 

distribution  over  the  wing.  The  lift  is  given  by 

N  N 

^  GT^On  tmmJ  %J  Qr%ork 


n=l 


span 


n=  1 


span 


where  F  is  the  circulation  distribution  at  each  section  of  the  wing.  The  lift 


n 


distribution  is  continuous  at  the  boundaries  of  the  jet.  Hence  at  the  points 
z  =  -h.  y  =  c^.c^....,C2j 


r  ,  ,(c,  ,-ih-0)  =  pr  ,  ,(c,  ,-ih+O) 

n,  Zv-r  lv~\  ^n,Z\>-Z  2w-l 


(77) 


r  ,  ,(c,  -ih+O)  =  ijT  ,  (c,  -i 

n,  Zw-r  Zv  '  n,  Zv  Zv 


ih-0) 


where  the  second  subscript  on  T  denotes  the  region  in  which  F  is  defined. 

The  circulation  distribution  F  (y)  at  each  section  is  defined  in  each 


n 


region  R  ,  for  v  =  0,1,...,  ZJ,  by 


rjyl  =  4bV.r„  „(y) 


(78) 


where 


A  -1 
o 


=  X  ^'zx+l.n"" 

A  .1 

V 

F  (y)  =  F  ,(y)  +A^'^^  +  ^  sin  { 2\.+l)4< 

n,  V  •'  n,v -1  •'  n,  o  ^  ZX.+l,n  v,n 


\=0 
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Here  the  dependence  of  the  angles  v|i  on  the  wing  section  chosen  has  been 

noted.  The  region  is  also  dependent  on  n,  being  defined  by 

c  >  V  ^  c  ,  where  c  =  a  -  Jl  -  h^  .  An  alternative  writing 
v.n  “  '  v+l,n  v,n  v  ^  n 


c  >  V  ^  c  ,  where  c 
V ,  n  ^  V  +1 ,  n 

of  78  is 


r  (y) 

n.  w  ’ 


ZX+1,  n 


sin  (  2X+1)4j  f  for  V  =  0, 1,  .  .  .  ,  2J 


with 


=  0 

o,  n 


Only  odd  multiples  of  the  angles  4<  are  taken  in  78,  which  relies  on  the 

assumption  that  the  geometric  incidence  of  the  wing  and  the  rotation  of  the 

slipstream  be  symmetric  about  the  mid-span. 

The  boundary  conditions  77  may  be  used  to  express  the  coefficients 

in  terms  of  the  unknown  coefficients  in  the  form 

o,n  i:A.+i,n 


(p-1)  1  V  -1 


>>  =  ,  y  A'='.y 


sin  for 

2\+l,  n  s  ,  n 


V  =  1,  3 ,  .  .  .  ,  2J-1 


v  =  2,  4,  .  . .  ,  2J 


8pbV 


The  total  lift  may  now  be  evaluated  as 


n= 1  V  =0 


v=»0  s=  0 


c  )A^®^ 

2v  -l,n'  o  ,n 


A  -1 
s 


V  .(8) 

'v  ,n  ^  ■'^2X+l,n\  ^ 

X=  0 

sin  2(X+1)^^^^  +  ^)  -  sin  2(X+1)3<^^^'  \1  j 

- -4(V-fT5 - - /J  J 
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C.  VELOCITY  POTENTIAL 


The  velocity  potential  of  a  vortex  element  at  an  arbitrary  height  with 
respect  to  a  single  jet  slipstream  has  been  determined  in  Section  II.  The 
method  used  there  involved  the  construction  of  the  image  system  of  the  vortex 
element  in  the  jet.  The  two-dimensional  part  of  the  solution  was  found  directly 
this  way.  A  similar  image  system  was  used  in  the  three-dimensional  part  of 
the  solution,  leaving  a  remainder  term  which  was  evaluated  analytically. 

The  velocity  potential  due  to  a  vortex  element  placed  in  a  system  of  multiple 
jets  will  be  constructed  using  this  image  method  as  follows.  The  vortex 


element  has  images  due  to  each  jet  separately.  Each  image  then  has  further 
images  due  to  the  rest  of  the  jets.  An  image  system  of  infinite  order  is 
thereby  generated.  An  approximate  solution  is  presented  which  takes  into 
account  the  primary  image  system  and  the  first  order  corrections  due  to  these 
images,  i.  e.  the  image-image  contributions.  However,  only  the  primary 
image  system  is  used  in  the  remainder  term  of  the  three-dimensional  part  of 
the  solution.  It  is  anticipated  that  this  term  is  already  small  2md  the  first 
order  corrections  are  neglected. 


The  two-dimensional  part  of  the  velocity  potential,  denoting  the  potential 


outside  all  the  jets  by  (p^  and  the  potential  inside  the  jet  by 

made  up  from  the  following  contributions.  The  factor  ^  r(T|)dTi 


is  omitted. 


(I)  For  a  vortex  element  r(T^)dT|  located  at  y  =  t]  ,  z  =  -h,  taken  outside  the 


Fj^(y-Ti,  z  +  h)  plus  a  sum  over  v  of  the  image  in  each  jet,  i.  e. 

*1  r  /  -  .  h 

[(Ti-a^)^+h^l^  L  (Ti-a^)^4h^  (Ti-a^) 

+  2(ti -a^  )hF2(  same  argument  as 


(Ti-a^l'+h' 


+  2(ti -a^  )hF2(  same  argument  as  F^ 
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This  latter  function  will  be  denoted  by  X  . 

^  V 

(l-€  ^)F,(y-  T^,  z+h)  plus  the  images  of  this  image  due  to  the  other  jets. 


pius  t 

i.  e.  a  sum  over  all  f  ^ 


\  of  (l-«  2)  • 


(II)  For  the  vortex  inside  the  jet,  f  ^  X. 


I 


(l-€  2)Fj^(y-Ti,  z -h)  plus  this  image  reflected  in  the  other  jets,  i.  e.  a 


sxim  over  all  v  /  (  of  (l-€  X^ 


(III)  For  the  vortex  inside  the  jet 


<P. 


i  ' 

<Pj.X 

^  all 


same  system  as  (II) 

F,(y-Ti,z+h)  and  the  single  image  -X.  .  The  vortex  element  reflected 
in  all  the  other  jets  gives  the  contribution  (1-t  ^)  X  summed  over 


The  construction  of  an  image  system  to  represent  part  of  the  three- 
dimensional  component  of  the  solution  is  accomplished  as  above  for  the  two- 
dimensional  part;  the  remainder  term  is  approximated  by  the  primary  image 
system  only.  The  complete  velocity  potential  may  now  be  written  down. 
Using  the  substitutions 


v-a  =  r  cos  8 

V  V  V 

v  +  a  =  r  cos  0 

/  v  -V  -V 

and  the  symmetry  about  the  mid-span. 


z  =  r  sin  0 

V  V 


z  =  r  sin  0  , 

-V  -V 


the  velocity  potential. 


outside  the 
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jets  is  given  by 


J  c. 


4Tr</>  =  (  r(ii)F(x,  y-n,  z+h)dTi  -  €  r(Ti)  [f(x,  y-ii,  z+h)+F(x,  y+ri,  z+h)]  dii 

°  -b  M  ''  c^. 


1=1  21 


^  J  . 

^  r(Ti)  ^  [x'(ri  )  +  X*  (-n  , )]  j  dri 


J  c 

-u-.,.y  f 

4.W  V 
1  =  1  21 


r(ii)[x  (n.a^)  +x  (-Ti,a^)]dr| 


J  J  c 


(-T1,  -a^)  +X  (-r;,  a^)  +  x'(ti,  -a^)]  dr| 


1=1  v=l  21 


J  GO 

y  r 

L  L  -jQ 

l=lm=0  ^ 


sin  kx 


X  sin  me^  K^(kq^)^  r  +  f  q^d  n 

^  ''‘^21  +  1  ''  ^2l  -1  ' 


P  21-1  , 

+  \  r(Ti)T'°'(k 


,’1  )dnj 


+  cos  m0 


21+1  ^Zl-\  ' 


r(Ti)C^^^(k,  Tl)dTl 


+  r  r(Tl)D^^^k.Tl)dTl 

^  r' 


"  ‘)r(W°Nk,.,)dn 

^  C-jj  ' 


+  sin  m6 


-21+2 
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+  cos  m0 


p^-2f 

r(Ti)Tj^’(k, 

‘^-2f +  1 

^-2f +  1 

+ 

"-2/ +2 

C  ' 

-2i 

^‘^-2f 

- 

•-2i  +  l 


where  X  (ti  t  a„)  is  the  extension  to  three-dimensions  of  the  function  x^ 
used  above,  that  is 


X{-n,a)=  - - 5- 

^  [(n-a^)Wl 


- 2 — T  - 2 — : 

(Ti-a^)Sh'^  (Ti-a^)Sh' 


h 

- 

-a  )  +h 


+  .!( Ti -a^)hG{ same  argument  as  F) 


The  velocity  potential,  <p .  .  ,  inside  the  X  jet  is  given  by 


□ 

^%,X=  I.  r(Ti)F(x,  y-r|,  z  +  h)dTi 


-b 

J  c. 


”  *  Z^^”*  2^ 


yp  ^  /  -i 

J  r(Ti  )[f(x,  y-Ti,  z+h)  +  F(x,  y+Ti,  z+h)]  dri 


i=l  2f 

J  ^  J 


+  (1-2) 


J  ^  ("n.  »v)  +  *  (ti  . -a^,)+X  (-T1,  a^,)+ X  (-Ti,-a^,)| 
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=>  I  tsj 
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D.  DOWNWASH  AND  BOUNDARY  CONDITIONS 


The  downwash  distribution  corresponding  to  the  velocity  potential 
for  the  multiple  jet  case  obtained  in  the  previous  section  may  be  obtained  in 
similar  fashion  to  the  simple  jet  case  (see  71).  The  integrals  involved  will 
be  those  listed  in  the  Appendix.  It  is  not  intended  to  write  down  the  downwash 
here  as  it  is  a  complicated  expression.  Assumption  of  a  slipstream  rotation 
of  the  type  described  in  III-D  will  lead  to  similar  boundary  conditions  to  73 
in  each  region  (see  also  Wu),  and  the  determination  of  the  coefficients  in 
the  expansion  79  for  the  circulation  will  proceed  as  before  from  the  solution 
of  a  system  of  linear  equations. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


Prediction  of  the  apanwiae  lift  distribution  on  a  wing  at  an  arbitrary 
height  and  high  angle  of  attack  extending  through  mxiltiple  jets  (propeller 
slipstreams)  has  been  accomplished  by  the  present  theory. 

Implementation  of  the  theory  requires  both  an  extensive  machine 
programming  effort  and  experimental  confirmation.  The  mathematical 
complexity  encountered  in  the  theory  requires  close  coordination  between 
the  analytical  and  machine  programming  efforts. 

It  is  therefore  recommended  that  the  following  two  (2)  steps  be 
undertaken: 

1.  Extensive  machine  calculations  -  a  programming  effort 
directed  towards  systematically  developing  a  generalized  calculation 
program  employing  the  wing  and  slipstream  geometry  as  inputs  should  be 
undertaken.  This  program  should  then  be  used  to  calculate  spanwise  lift 
distributions  of  a  broad  spectrum  of  arrangements  to  (1)  delineate  significant 
trends  due  to  non-uniform  flow  effects  and  (2)  provide  results  for  comparison 
with  experimental  data. 

2.  Systematic  Experimental  Program  -  wind  tunnel  tests  to 
confirm  the  analytical  results  and  determine  real  fluid  effects  should  be 
conducted.  The  tests  should  be  designed  to  yield  accurate  spanwise  lift 
distributions  of  wing- slipstream  arrangements  selected  from  the  analytical 
study  above  to  yield  a  meaningful  confirmation  of  the  theory. 
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APPENDIX 


A  series  of  integral  functions  denoted  by  were  introduced  in 

Section  III-  C,  These  integrals  will  be  defined  in  the  following  and  the  param¬ 
eters  used  in  their  definitions  will  be  noted  as  arguments.  A  number  of 
supplementary  functions  are  used  which  are  given  by 

(y  -  b  cos  9)^-  (H  -|-  h  )  ^ 

A(0)  =  - X - - - sin  0 

[(y  -  b  cos  9)  +  ^ 
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The  integrals  are  defined  by 

®  n  a 

,1^  (y.b.H  .h  ,\)  =  r^A(e)  sin  Xe  de 
1  O'  ^  r  n  J 

a 

^1^  (y.b.H^.h^.X^.x^.X)  =  (X^-  x^)j'^A(e)C(0)  sin  Xe  dG 


r 


3 


3!'^  {y.b,H^.h^.X^.x^.X)=  (X^-x^)(H^+h^)‘'\  B(e)C"(e)  sin  xe  dG 

^  (y.b.H^.h^)  =  r^A(e)de 
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13^  =  Iq  Dj(e)Ej(e)E2(e)  de 
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